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Understanding the Optical Potential in HBT Interferometry 



I. INTRODUCTION 
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The validity of using a pion optical potential to incorporate the effects of final state interactions 
on HBT interferometry is investigated. We find that if the optical potential is real, the standard 
formalism is modified as previously described in the literature. However, if the optical potential is 
complex, a new term involving pion emission from eliminated states must be included. The size of 
such effects in previous work by Cramer and Miller is assessed. 

PACS numbers: 25.75.-q, 25.80.Ls, 13.85.Hd 
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C/^ . The space-time structure of the "fireball" produced in the collision between two relativistically 
CN " moving heavy ions can be investigated by measuring two-particle momentum correlations between 
^ pairs of identical bosons. The Bose-Einstein enhancement of the coincidence rate at small momentum 

^ differences depends on the space-time extent of the particle source. This method of investigation, 
' I . called HBT interferometry, has been applied extensively in recent experiments at the Relativistic 
g Heavy Ion Collider (RHIC) by the STAR and PHENIX collaborations. See the reviews[l-4]. 

The invariant ratio of the cross section for the production of two pions of momenta pi, p2 to the 
' — ' product of single particle production cross sections is analyzed as the correlation function C(pi, P2). 
^ . We define q=Pi-p2 and K=(pi+p2)/2, with as the component perpendicular to the beam 
^ [ direction. (We focus on mid-rapidity data, where K = K^.) The correlation function can be 
^ parameterized for small q as C(q, K) — 1 ^ A exp (— -Rq^'o ~ ^sls ~ ^lQl) ~ ~ ^o^o ~ ^sls ~ 
10 RlIl) {liRi ^ 1); where O, S, L represent directions parallel to Kt, perpendicular to both Kt and 
. the beam direction, and parallel to the beam direction [5]. Early [6] and recent [3] hydrodynamic 
^ I calculations predicted that a fireball evolving through a quark-gluon-hadronic phase transition would 
QQ ' emit pions over a long time period, causing a large ratio Rq/Rs- The puzzling experimental result 
O that Ro /Rs ~ 1 [7] is part of what has been called "the RHIC HBT puzzle" [8] . Another part of the 
>■ \ puzzle is that the measured radii depend strongly on the average momentum K, typically decreasing 
[ in size by about 50% over the measured range, showing that the radii are not simply a property of 
5^ ' a static source. The medium at RHIC seems to be a very high density, strongly interacting plasma 
[9], so that any pions made in its interior could be expected to interact strongly before emerging. 
Thus one expects that the influence of the interactions between the pion probe and the medium, as 
well as flow and other effects, must be taken into account when extracting the radii. 

We studied the effects of including the pionic interactions in previous work [10, 11]. Distorted 
waves, instead of plane waves, were used to represent the pion wave functions. The resulting for- 
malism is called the Distorted Wave Emission Function (DWEF) formalism because the emission 
function used to describe the space-time extent of the emitting system is dressed by the final state 
interactions. We found that it is possible to simultaneously describe the measured HBT radii and 
pionic spectra by including the effects of pion-medium final state interactions obtained by solving 
the relevant relativistic wave equation. These interactions are so strongly attractive that the pions 
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act as essentially massless objects inside the medium. The medium acts as if it is free of the chiral 
condensate that is the source of the pion mass, and therefore acts as a system with a restored chiral 
symmetry. Other solutions of the HBT puzzle have been proposed. See the review [4] . 

Some theorists have questioned whether waves, produced by incoherent sources, unaffected by 
final state interactions, interfere with those that are affected by final state interactions. That this 
interference occurs was demonstrated at least as early as 1979 [12], again in the nineties [13-16], 
and there have been two recent publications confirming that conclusion [17], [18]. The former derive 
a general formula for the correlation function of two identical particles including multiple elastic 
scatterings in the medium in which the two particles are produced. Numerical results for the case 
of soft final state interactions are presented. Ref. [18] includes the effects occurring when emitted 
particles undergo multiple scattering with medium particles. Using the Glauber theory of multiple 
scattering at high energies and the optical model at intermediate energies, it is found that multiple 
scattering leads to an absorption. 

Despite this progress, several conceptual issues remain. These include understanding: the meaning 
of the imaginary part of the optical potential, the role of the energy dependence of the optical 
potential [19], and the relationship between the sources producing the pions and the optical potential. 
Therefore, we find it worthwhile to re-investigate the effects of quantum mechanical treatments of 
final state interactions. 

Our procedure is to repeat the derivation of Ref. [5] using a simple Lagrangian. First, the original 
plane wave treatment is reproduced using our notation. Sect. II. Then the effects of a real opti- 
cal potential are incorporated. Sect. III. The result is a re-derivation of of the DWEF formalism. 
However, a deeper understanding is needed to correctly account for the effects of a complex optical 
potential. This can only be incorporated using a coupled channels formalism. Sect. IV. We find 
that including the complete effects of an imaginary optical potential requires a modification to the 
DWEF formalism that is presently incalculable. However, the optical potential used in [10, 11] was 
dominated by its real part. In particular, in Sect. V we find that setting the imaginary part optical 
potential to zero does not significantly change our description of the data. Sect. VI is reserved for a 
summary and discussion. 



The space-time extent of a source of pions can be inferred by measuring the pionic correlations 
known as the Hanbury Brown-Twiss effect [20, 21]. The correlation function function C(p, q) is 
defined to be 



where P{pi, ■ ■ ■ Pn) is the probability of observing pions of momentum pj all in the same event. The 
identical nature of all pions of the same charge cause C(p, p) = 2. The width of the correlation 
function is related to the space-time extent of the source. 
A state created by a random pion source |7^) is described by [12] 



factor that takes the chaotic nature of the source into account, and c^(p) is the creation operator 
for a pion of momentum p.. In particular, an average over collision events gives 



(7*(t)7(t')) = S{t - t'), {Y{hh*{t2h{t3h{U)) = S{h - h)S{t2 - U) + 5{h - U)S{t2 - h). (3) 



II. THE PRATT FORMALISM 




(1) 




We note that as written, the state \ri) is not normahzed to one. However, the normahzation constant 
will divide out of the numerator and denominator of the correlation function. Therefore we do not 
make the normalization factor explicit here, but note that it enters when we calculate the pion 
spectrum. 

For ■0 and its time derivative to obey the Heisenberg commutation relation one has 

v/^[c(p),ct(p')]=5(^)(p-p'). (4) 

Furthermore, we define 

77(p,i) = j (fxe-'P-''r]{x). (5) 
The state \f]) is an eigenstate of the destruction operator in the Schroedinger representation, c(p): 

cip)\r)) = I dte'^^^'^mri)- (6) 

The correlation function is 

^/ X ^ (?7|ct(p)ct(q)c(q)c(p)|?7) 
^P'^^ (r?|ct(p)c(p)|r^)(r^|ct(q)c(q)|r^)- ^ ^ 

The use of Eq. (3) and Eq. (6) in the numerator of Eq. (7) yields 

(^|ct(p)ct(q)c(q)c(p)|r^) = {'n\c\p)c{pM{v\c\(i)c{(iM + |(r/|ct(p)c(q)|r;)|2. (8) 
Furthermore 

(r^|ct(p)c(q)|r^) = / dteM-i{E, - E,)t]'^^^^^^^^^^^ . (9) 
The quantity g{x,p) is denoted the emission function and is defined as 

g{x, p)=J d'x'r,*{^ + ^x', t)r,{^ - ^x', t)e*P-', (10) 

so that 

"^'^ ^((y + y')/2,t,p)e-^P-(^-^') = r^\y,tUy\t). (12) 



/ 



(27r) 



The second expression appears in the right-hand-side of Eq. (9) (if one uses Eq. (5)) so that we may 
write 

I \ \( \ ( w \ /"^4 exp[-i(p - g) • x] (p + g). 

(77|ci^(p)c(q)|r7) = y d^'x ^-^ ^(x, — ). (13) 

Using Eq. (13) with p = q shows that the function g{x^ p)/Ep is the probability of emitting a pion 
of momentum p from a space-time point x. Using Eq. (8) and Eq. (13) in Eq. (7) gives the desired 
expression: 

, _ J d^x d^x'g{x,^K)g{x\^K)exp[-ik ■ {x - x')] 
0(p,qj-i+ Jd^xd^x'g{x,p)g{x',q) ' ^ ' 
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where K = p + q and k = [Ep — Eg, p — q), and the factors of have canceled out. 

From a formal point of view, a key step in the algebra is the relation between the Heisenberg 
representation pion creation operator ^^(x) and its momentum-space Schroedinger representation 
counterpart c^{p) that appears in Eq. (2): 

^ I d'p cHp) e^^^' (15) 

i^ix)^ I d'p c{p) 6-^"=^' (16) 

The operators c^i^p) (c(p)) are coefficients of a plane wave expansion for ip^i^x) (ip^x)), with the plane 
wave functions e*P'^/(27r)^/^ being the complete set of basis functions. However, one could re-write 
ip^{x) (^(a;)) as an expansion using any set of complete wave functions. We shall exploit this feature 
below. 



III. DISTORTED WAVES - REAL POTENTIAL 

We represent the random classical source emitting pions that interact with a real, time- independent 
external potential U by the Lagrangian density: 

-C^i:\-d'^ + U + m'^)il: + j{x)il^. (17) 

The current operator j{x) is closely related to the emission function g, [12]. In this Lagrangian the 
terms U and j{x) are independent. Thus the relation between the emission function and U derived 
in [22] need not be satisfied. 
The field operator ij)^ can be expanded in the mode functions ■i/^p"'' that satisfy: 

(_V2+Z^)^(-)(x)=pV(-)(x). (18) 

These wave functions obey the usual completeness and orthogonality relations 

J d3pV'(-)*(x)V^(-)(y) = 5(3) (x - y) (19) 

J d3xV^(-)*(x)^(7)(x) = S^'\p - pO, (20) 

so that one may use the field expansion 

i^ix) = J d3pV^(-)(x,t)e-^-*d(p), (21) 

with S{p) being the creation operator for pions of momentum p in the basis of Eq. (18). The 
expansion Eq. (21) assumes that U produces no bound states. If so, one the integral term would be 
augmented by a term involving a sum over discrete states. 

The availability of mode expansions when distortion cfi^ects are included means that the simplifi- 
cation of the correlation function can proceed as in the previous section. We again use Eq. (2) and 
Eq. (3). The use of the field expansion Eq. (21) enables a generahzation of the function ri{x): 



ri(x)^ / £p^|;^^>{^,t)fj(p,t), 



(22) 



with 

fl{p,t) = J rf3x^(-)(x,t)r?(a;), (23) 

so that 

\V) = cxp[ J d'pdtfj{p,tMt)d\pm. (24) 

The abihty to obtain a relation between the fj{p,t) and 'r){x) rests on the relations Eq. (19) and 
Eq. (20). 

The state \ri) is an eigenstate of d{p). Thus the result 

rrn = 1 + \{v\dKp)d{q)\v)\' 
^P'"*^ ^ {rj\d^{p)d{pM{v\d^{ci)d{ci)\vy ^ ' 

very similar to Eq. (7), is obtained. We need the matrix elements appearing in the numerator and 
find 

{r,\d^{p)dici)\r,) = J d'xd'x' ^^^^ i^'HxH^', t) ■ (26) 
and the use of Eq. (12) allows us to obtain 

{V\d\p)d{q)\r}) = ^J citci3a;dV|^e^*(^^-^-)e-P'-V(-)(x + x72)^(-)*(x-x^ 

(27) 

This result, which can be applied for p q and for p = q, specifies the evaluation of the correlation 
function of Eq. (25) with the result 

S(K.k) , , 

where 

S{K, k) = j ciWx'^^e'*(^^-^p')e-'P'-"V(,-)(x + x72)V^^-)*(x - x72)^(x, p'), (29) 

and 

S{p) ^ J ciWa;'^e-P'-V(-)(x + x72)V'i-)*(x-x72)5(x,p'). (30) 

This expression is also the one that appears in the DWEF formalism [10, 11]. One could use either 
Eq. (14) or Eq. (28) to analyze data, but the extracted space time properties of the source rj{x) 
would be different. 

We need to comment on the possible momentum and energy dependence of the optical potential. 
The completeness and orthogonality relations are obtained with any Hermitian U which can therefore 
be momentum dependent, but not energy dependent. As explained in Sect. 5 (Eq. (43)) of Ref. [11], 
the real part of the potential can and should be thought of as a momentum-dependent, but energy- 
independent potential. If there were true energy dependence a factor depending on the derivative 
of the potential with respect to energy, [19], would enter into the orthogonality and completeness 
relations. 
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IV. COUPLED CHANNELS 

The optical potential used in previous work [10, 11] is complex. Using the necessary completeness 
and orthogonality relations to relate r]{x) to ?7(p, t) requires the use of a real potential. Therefore 
one needs to investigate possible corrections. 

The optical potential or pion self-energy is an effective interaction between the pion and the 
medium. The medium is not an eigenstate of the Hamiltonian, but rather of Hq, which is the full 
Hamiltonian minus the Hermitian operator representing the pionic final state interactions. Elimi- 
nating the infinite number of possible states of Hq and representing these by a single state leads to 
a self-energy that is necessarily complex. Our procedure here is to specifically consider the infinite 
number of states of the medium, obtain a Lagrangian density that involves Hermitian interactions, 
and derive the optical potential formalism and any corrections to it. 

Let Pn denote a projection operator for the medium to be in a given eigenstate of Hq, n. These 
obey 



Pn 1) PnPm ^n,mPn- (^■^) 



For the case of 7r-nuclear scattering, n would represent the nuclear eigenstates. Here n represents 
states of the medium in the absence of its interactions with pions. The correlation function is now 
given by 

where i-*n(p) is the probability for emission of a pion of momentum p from the medium in a state 
n. Similarly -P„(p, q) is the probability for emission of a pair of pions of momentum p, q from the 
medium in a state n. The sums over n account for the inclusive nature of the process of interest. 
It is convenient to define the product of the field operator with the projection operator P„: 

1pn{x) = 1p{x) Pn, (33) 

with 

^{x) = Y,^n{x), (34) 

n 

using the complete nature of the set n. The Lagrangian density is given by 

-£ = ^ • 94 + ^ {{ml + Ml)5nm + Unm) V^m + J] 3n{x)i^n{x), (35) 

n n,m n 

where 

Unm^Kn=i^<)nm (36) 

and U is the Hermitian interaction operator and M^, the m matrix element of the diagonal operator 
M^, represents the effects of the different energies of the states labeled by m. The field operator 
can be expanded in the mode functions V'pT^- 

J2 Unm{^)^U{^: t) = {p' + - - W„„(x))4;) (x, t). (37) 



Here the potential lA is taken as a local operator in the position space of the outgoing pion. 

To see the correspondence between the formulation of Eq. (35) and Eq. (37), let ipi correspond to 
the field operator (and state) of the previous section and solve formally for ip^X terms of ip'^ {*. It 
is convenient to define the operator U with matrix elements given by 

Un,n' = (1 — ^n,l)(l " ^n',l)^n,n' (38) 

Then 

where (V^ — M^)„m oc 5n,m, and is an operator giving when acting on the state n. Then 
rewrite Eq. (37) in terms of ip^J as 



(40) 



The complex object 



a non-local operator in coordinate space, can be identified with the optical potential, given by the 
operator V{Z) as a function of a complex variable Z: 

m,n^l \ ^ ' ' / m,n 

We proceed by employing Eq. (34) and Eq. (35) to compute the correlation function. The solutions 
of Eq. (37) form a complete orthogonal set: 

E / ^'p'^^'-^i* (^X:^ (y) = '^^'^ - y) (42) 

E / '^'^ V'i^n^* WV'i'iW = ^^'Hp - P')- (43) 

The field expansion is now 

^(x) ^jd'pJ2 a(p)P„^(;)(x)e-^^-*, (44) 



so that 



|r;) = exp[5] y d^a;r;,(a;)7(t) / a\p)Pn^^J*{^)e''^-'] J] |0, m), (45) 

n m 

where the state |0,m) is the pionic vacuum if the medium is in the state m, and rin{x) represents 
the source for the state n. These state vectors obey the relations 

(0,n|0,m) = 6n,m = (0, n|P„|0, m). (46) 
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Define 

Vn{p,t) = J d^X7]n{^,t)ijl->{^), (47) 

SO that 

\7j) = exp[ J d'pdt^{t)^7]n{p,t)a\p)Pne'''^'] |0,m), (48) 

n m 

a{pM = fdt 7(t) Yl ^^^Pne'^'^'^lv)- (49) 
The emission probabihty is given by 

EpE, J2iv\aKp)Pna{ciM = E / dtd'xd'yrjli^, t)r)^{y, t)V'i-)*(x)V'i-)(y) e^^^^"^")* (50) 

n n 

or using Eq. (12) 

EpE,{rj\a\p)a{<i)\v) = J2 f ^^^'^^'^ / ^n((x + y)/2, t, p')e-P'-(--^)V^(;)*(x)V^(-)(y) e^(^''-^^)*, 



(51) 



where 



^„(x, p) = y" d^x'<(x + ix', i)77n(x - ^x', t)e^P '''. (52) 

If picnic final state interactions are ignored, the term gf„ enters and this may be identified with 
the emission function, g of previous sections. 

The expression Eq. (51) is the same as Eq. (27) except that now we sum over the channels n. These 
sums may be expressed in terms of the optical model wave functions of Eq. (39). The term of Eq. (51) 
with n = 1 corresponds to the DWEF formalism, and the terms with n > 1 are corrections. We 
provide an example of a correction term. Suppose part of the imaginary part of the optical potential 
arises from a pion-nucleon interaction that makes an intermediate A. Then a term corresponding to 
one of n > 1 involves the emission of a pion from a nucleon that makes an intermediate A. 

It is difficult to assess the importance of the second term in a general way. The only obvious limit 
is that if states with n > 1 are not excited then Im{V) of Eq. (41) must vanish. Conversely, if 
Im{V)=0, the states n > 1 must be above the threshold energy and the propagators that appear in 
the correction terms correspond to virtual propagation over a small distance with limited effect. 

V. NUMERICAL ASSESSMENT OF THE EFFECTS OF Im{V) IN REFS. [10, 11] 

We proceed by assessing the possible importance of the correction term for the work of [10, 11] 
by seeing what happens if the optical potential is taken to be purely real with no imaginary optical 
potential. A variety DWEF fits are performed, see Table I. In [11] the imaginary part of the optical 
potential as represented by the term W2 is about one tenth of the real potential. It is therefore 
possible that, in the limit that Im{w2) = 0, there would be no significant correction term, so we try 
to understand if removing the imaginary part of the optical potential can be done without degrading 
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TABLE I: Four parameters sets obtained with slightly different procedures. These parameters are defined in [11]. The values 
of represents the accuracy of the description of the data. 



T{MeV) rif 


Ar(/m/c) 


Rwstfm) 


awsifm) 




W2 To{fm/c) 


A?7 






121 1.05 





11.7 


1.11 


0.495 


0.762 +0.0001i 9.20 


70.7 


139.57 


300 


162 1.22 


1.55 


11.9 


1.13 


0.488 


1.19+0.13i 9.10 


1.68 


139.57 


117 


121 1.04 


1.5 


11.7 


0.905 


0.564 


0.595 +0.0001i 8.85 


70.7 


139.57 


451 


144 0.990 


2.07 


12.57 


0.876 


0.0001 


0.0001+ O.OOOli 6.85 


oo 


83.5 


1068 



the quality of the fit. The results are shown in Figs. 1 and 2. An example of the previous calculations 
[10, 11] is shown as the green dotted curve (second hne of Table I). The red sohd curve (first hne 
of Tabic I) shows the result of setting the imaginary potential to a vanishingly small value. This 
results in only a slightly worse description of the data. The changes in the imaginary part of the 
optical potential W2 are largely compensated by a reduction of the temperature from about 160 MeV 
to about 120 MeV. We also point out that the length of the flux tube as represented by Ar] is vastly 
increased, providing greater justification to our previous procedure of taking the length of the fiux 
tube to be infinitely long in the longitudinal direction. However, the emission duration is reduced 
to fm/c, which is similar to the results of the blast wave model [25]. This means that all of the 
pionic emission occurs at a single proper time. This value justifies the use of a time-independent 
optical potential, but does seem to be difficult to understand because some spread of emission times 
is expected for a long-lived plasma. The results shown by the blue dashed curves (third line of Table 
1) arc obtained with fixing the emission duration to 1.5 fm/c, which is our previous value [10, 11]. 
The description of the spectrum is basically unchanged but the radii are less precisely described. 
The violet long-dashed curves (fourth line of Table I) show the DWEF fit using a vanishing optical 
potential. This does not give a good description of the momentum dependence of the radii and is 
associate with the largest deviation between our calculations and the data as represented by the 
values of Table 1. 

It is clear that the precision of our description of the data is improved by including the imaginary 
part of the optical potential. However, this is a quantitatively but not a qualitatively important 
effect. It is also true that including the real part of the optical potential is a qualitatively important 
effect. These results suggest that the correction terms embodied by the terms with n ^ 1 oi Eq. (51) 
are not very important, but non- negligible. It is also possible that an optical potential with a 
different geometry than the volume form that we have assumed might be able to account for the 
the neglected terms. However, an accurate assessment would require the development a theory that 
involves dealing with explicit models for gn,jn and U. 

VI. SUMMARY AND DISCUSSION 

It seems clear from previous work including [12]- [15] and the present Sect. HI that final state 
interaction effects on HBT interferometry are appropriately included by solving quantum mechanical 
wave equations. However, if the optical potential has an imaginary part, there is an additional effect, 
embodied in Eq. (51) that needs to be included when computing the emission probabilities and 
correlation function. Thus the effects of strong quantum opacity must be accompanied by additional 
pion emission from the states eliminated in the construction of the complex optical potential. In the 
work of [10, 11] the real part of the optical potential is very important and the imaginary part of 
the optical potential is a small effect. However, obtaining a similarly accurate reproduction of the 
pionic spectra and HBT radii without this imaginary part causes the emission temperature to drop 
from about 160 MeV to 120 MeV and the fitted emission duration time to drop to 0. This indicates 
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FIG. 1; Color online. Computed pionic spectrum. Red upright triangles: tt" spectrum (STAR) Green inverted triangles; tv^ 
spectrum points (STAR) [24] Red solid line: DWEF fit with vanishing imaginary part of the optical potential, first line of 
Table I. Green dotted line: DWEF fit including search on the imaginary part of the optical potential, second line of Table 
I. Blue dashed line(almost entirely covered by the red solid curve): DWEF fit with vanishing imaginary part of the optical 
potential,Ar = 1.5 fm/c, third line of Table I. Violet long dashed line: DWEF fit including search on Hn, setting the optical 
potential to essentially 0, fourth line of Table I. 



that, in our model, either the final state interactions occur in the later times of the collision and that 
the emission occurs at only one proper time, or the inclusion of emission from the states eliminated 
in the construction of the complex optical potential is necessary. 
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